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Abstract 

We discuss the cosmological evolution of a 3-brane Universe in the presence of energy 
influx from the bulk. We show that this influx can lead to accelerated expansion on the 
brane, depending on the equations of state of the bulk and brane matter. The absorption 
of non-relativistic bulk matter by the brane at an increasing rate leads to a small positive 
acceleration parameter during the era of matter domination on the brane. On the other 
hand, the brane expansion remains decelerating during radiation domination. 


Introduction: Recent astronomical observations indicate that the expansion of our Universe is 
accelerating [1], The physical mechanism that drives this cosmological acceleration has not been 
established yet. The most popular explanation relies on the presence of a small cosmological 
constant or dark energy. However, in order not to disturb successful ingredients of standard 
Big Band cosmology such as nucleosynthesis, the contribution of the dark energy to the total 
energy density must be strongly constrained. As a result, the dark energy can be significant 
only during recent times. A convincing explanation of this “cosmic coincidence” has not been 
given yet. Most scenaria identify the dark energy with the vacuum energy of a time-dependent 
field with carefully chosen dynamics [ 2 ]. 

In this letter we would like to explore another possibility, that does not require the presence 
of a field. The mechanism can be demonstrated by considering the equation 

p + ^H{p + p)=T, (1) 

that describes energy conservation in an expanding Universe with Hubble parameter H, in the 
presence of a source term T. The standard conservation equation with T = 0 has stationary 
solutions with p = 0 only if p + p = 0. This implies the presence of vacuum energy with 
equation of state p = —p, that leads to the standard inflationary scenario with constant H and 
an exponentially increasing scale factor a. 

However, if T 7 ^ 0 stationary solutions are possible for a more general class of equations of 
state. For example, one could have constant H, p, p, such that 2>H {p + p) = T. The resulting 
cosmological solution has an exponentially increasing scale factor and positive acceleration. It 
has been termed “steady-state Universe” [3]. Interesting behaviour, with positive acceleration, 
can be obtained even if T is not constant, but has some dependence on the energy density p or 
the scale factor a. 

In the framework of a four-dimensional world, the energy source term T does not have a 
simple origin. However, if more than four dimensions exist, our observable Universe may be 
identified with a 3-brane embedded in the higher-dimensional world. In this case, eq. (1) would 
refer to the matter localized on the brane. The source term T could have the simple interpreta¬ 
tion of energy density falling onto the brane from the extra dimensions. The energy-momentum 
tensor in the higher-dimensional world is conserved and no source terms exist. However, an 
effective source term is generated for the matter localized on the brane. In the following we 
shall discuss how such a scenario can be realized. 

The framework: We work in the framework of the Randall-Sundrum (RS) model [4]. It is 
described by the action 

5 = J d^x ^ (m^R - a + + I (-U + £r*) , ( 2 ) 

where R is the curvature scalar of the five-dimensional metric qab, A, H = 0,1, 2,3,4, and 
with a, f3 = 0, 1, 2 ,3, is the induced metric on the 3-brane. A is the bulk cosmological constant, 
while the quantity V includes the brane tension as well as quantum contributions to the four¬ 
dimensional cosmological constant. We consider the commonly used ansatz for the metric [5, 6 ] 

z)dt^ -[- a^{t, z)'yijdx^dx^ + b^{t, z)dz^, (3) 

where 'jij is a maximally symmetric 3-dimensional metric. We use k = —1,0,1 to parametrize 
the spatial curvature. 
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We decompose the energy-momentum tensor into vacuum and matter contributions in the 
bulk and on the brane 


rpA 

^ c 

= T^c +Ti +T^c +Ti ^ 

v,o m,b m^B 

(4) 

Tjj 

VjO 

= Mdiag(-I/,-E,-I/,-E,0) 

(5) 

> 

= diag(—A, —A, —A, —A, —A) 

(6) 

T^c , 

m,6 

S(z) 

= ^ diag( p,p,p,p,0), 

(7) 


where p and p are the energy density and pressure on the brane, respectively. The behaviour 
of is in general complicated in the presence of flows, but we do not have to specify it 

further in this work. 

We are interested in the Einstein equations at the location of the brane. We indicate by 
the subscript o the value of various quantities on the brane. It is convenient to work in a 
coordinate frame in which = no = 1. This can be achieved by using coordinates such that 
6 (t, z) = 1. This is always possible as can be checked by considering the “two-dimensional” part 
of the metric: z)dt^ + z)dz‘^. Through an appropriate coordinate transformation 

we can always set this part in the form z)dt^ + dz'^d We can then redefine the time 

coordinate so that no = 1 on the brane. We emphasize that our assumptions for the form of the 
energy-momentum tensor become now specific to this coordinate frame. In particular, the brane 
is identified with the hypersurface z = 0 in this frame. In this sense, our discussion becomes 
more restricted. The gained advantage is that significant progress can be made without solving 
explicitly the Einstein equations in the bulk in the presence of flows. 


The Einstein equations: Eor b{t, z) = no = 1 we find [5]-[8] 


p + 3^{p + p) = -2T% 

do 


- { — H—^ — 


12M3 


144M6 


{V{3p - p) + p{3p + p)) - 


where r 4 ,T^ are the 04 and 44 components of evaluated on the brane. 

We are interested in a model that reduces to the RS vacuum [4] in the absence of matter. 
In this case, the first term in the r.h.s. of eq. (9) vanishes. A new scale is defined through 
the relations V = —K/k^g = \2M^kgig. 

We can rewrite eqs. (8), (9) in the equivalent form [7, 8] 


p + 3Ho{p + p) = 


tt 2 _ “o 

2 

«o 


ic + AHoX = 


^ + AHo4> = - 


-2T\ 

1 

144M6 

1 

36M6 * 

1 


{p^ + 2Vp) - -^-kx + (/> + A 


ip + V)T\, 


- zHoTi, 

3M3 


^The conformally flat form c?{z,t){—dt^ ^dz^) is more commonly used in the literature. The important point 
is that, using the reparametrization invariances, a two-dimensional metric can be put in a form that depends only 
on one function of t,z. 
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where A = (A + y^/ 12 M^)/ 12 M^ is the effective cosmological constant on the brane. The 
functions X) defined through eqs. (12), (13). In the RS model A is zero, the value we shall 
use in the rest of the paper. We have assumed the orbifold symmetry 2 : — 2 , so that 2T\ is 

the discontinuity of the 04 component of the bulk energy-momentum tensor at the location of 
the brane. 

In the low-density region, in which p <C R, we may ignore the term ~ p in the above 
equations compared to V and define Mpi = 12M^/V = M^/kpg. The ratio of the terms in the 
r.h.s. of eqs. (13), (12) is of order 

n iikk. 41 ,/z, ( 14 ) 

We assume that p/V is sufficiently small for this ratio to be small at all times of interest^. Then 
4> can be set to zero and omitted from our considerations. As a result, eqs. ( 10 )-( 12 ) can be 
written as 


p-|-3(l + w)HoP 



X + 4.H0X 


1 

to 

(15) 

p k 

QMpi ^ al 

(16) 

CO 

(17) 


where p = wp. The cosmological evolution is determined by three initial parameters (p*, Uj, Xi: 
or alternatively pi, ai, ai), instead of the two {pi, a*) in conventional cosmology. The reason is 
that the generalized Friedmann eq. (11) (or (16)) is not a first integral of the Einstein equations 
because of the possible energy exchange between the brane and the bulk. In the above equations, 
we recover the “mirage” or “Weyl radiation” component X; found in studies of the cosmological 
evolution in the presence of energy outflow [ 8 , 9]. 


The bulk energy-momentum tensor: The determination of the cosmological evolution on 
the brane requires information on the form of T 4 . This can be determined exactly only through 
the solution of the Einstein equations in the bulk, a formidable task in the case of energy flows. 
However, the qualitative behaviour of T 4 can be inferred from the form of the energy-momentum 
tensor. For a perfect fluid with five-velocity it has the form 


rj-^AC 


— + (Pb + PbW^U^ ■ 


(18) 


For a fluid falling onto the brane with velocity along the fifth dimension, and within our 
ansatz for the metric, we have = (1 — W = 0, at the location of the 

brane. As there is no expansion along the fifth dimension, we expect that the main effect of the 
expansion along the remaining three spatial dimensions is to dilute the fluid energy density and 
reduce the pressure. These should fall with a certain power of the scale factor, determined by 
the equation of state of the fluid. 

^It is obvious that the omission of (j> is not justified when the energy influx stops and T 4 becomes zero. 
However, we assume that at this point T\ is also very small, so that <j) gives a negligible contribution to eq. (11). 
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In order to determine this power we may consider the conservation of the energy-momentum 
tensor in the absence of energy flows. For the 00 component of at the location of the 

brane we find 

n = -(-3n+n), (19) 

do ^ ^ 

where we have made use of 6o = Uo = 1- For a non-relativistic gas of bulk particles with 
= diag(—/ 3 n, 0, 0,0, 0), we obtain pp ~ a~^. For a relativistic gas with = 

diag(—pp, pp/4, pp/4, pp/4, pp/4), we obtain pp = 4pp ~ Oo . For a gas with zero pressure 
along the fifth dimension and = diag(—pp, pp/3, p^l^i Pb/^i 9); obtain pp = 3pp ~ 

in,S 

a“^. Finally, for the case of vacuum energy with = diag(— pp, — Pp, —Pr, —Pn^ ~Pb)^ 

we obtain pp = —pp ~ const. 

In analogy with the above, for a non-relativistic gas of bulk particles with constant we 
expect T 4 ~ a~^ at the location of the brane, while for an isotropic relativistic gas T 4 ~ ao . 
liT^ originates in the vacuum energy of a bulk scalar field we expect T 4 ~ const. Other types 
of behaviour are also possible in scenaria in which the velocity U 4 varies with time. 

We parametrize the dependence of T 4 on Oq as 



For energy influx, we have T 4 < 0, T > 0. We assume implicitly that the form of T 4 originates 
mainly in the bulk dynamics, determined through the bulk equations of motion. The presence 
of the brane generates only a small perturbation to the bulk evolution. This is not expected 
to be always the case. For example, if the energy density on the brane exceeds a certain value 
one would expect the influx to stop. This means that T 4 could vanish or even change sign 
at the location of the brane, even though it could remain unaffected far from the brane. The 
complicated dynamics associated with such phenomena is beyond the scope of this work. For our 
purposes we shall rely on the simple ansatz of eq. ( 20 ) that accounts mainly for the dilution of 
energy density through expansion. The variable T in eq. (20) is taken to be independent of the 
brane parameters, an approximation that is valid only for hnite ranges of the brane evolution. 

In this work we do not specify the mechanism that is responsible for the transfer of energy 
onto the brane. However, the energy influx seems a natural phenomenon in models in which the 
brane particles are identified with light modes of the five-dimensional theory that are localized 
on some defect [10]. In a fluctuating system, in which the light modes are not significantly 
populated, one would expect energy to be transferred from the massive modes to the light ones 
through interactions. As the light modes are localized on the defect, energy is expected to flow 
towards it. 

The substitution of eq. (20) into eqs. (15)“(17) results in a closed system of equations. In 
the following we neglect the spatial curvature and set /c = 0 in eq. (16). 

Solutions: We are interested in solutions that describe accelerating eras in the cosmological 
evolution of the brane, without the presence of an effective cosmological constant. For this 
reason, we concentrate on the range 0 < rc < 1/3 for the parameter that determines the 
equation of state of the brane matter. We also consider the range 0 < g < 4 for the parameter 
appearing in eq. ( 20 ), which is directly related to the equation of state of the bulk matter as 
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we explained above. In a stationary state, in which particles with constant velocity along 
the fifth dimension fall onto the brane, we expect q = 3, 15/4 or 4. Slow variations of ^4 can 
be modelled by allowing non-integer values of q. For example an increase of |r; 4 | for infalling 
massive matter should lead to an increasing parameter T in eq. (20) or, alternatively, to a value 
of q smaller than 3. We do not specify the details of the mechanism through which energy is 
transferred from the bulk to the brane. This can be done only within a scenario that includes 
a dynamical localization mechanism for the brane matter. 

For the range 0 < tc < 1/3, 0 < (7 < 4 the system of eqs. (15)-(17), (20) has solutions of the 
form 


p{ao) _ Cl 
QMpi 

x{ao) = 

^0 

( 21 ) 

s = 

2ql3 

( 22 ) 

Cl 

4 — s 

(23) 

C 2 

3(1 + w) — s 

{Ci + C2f^ = 

l-3w ^ 

(24) 

[3(1 + re) - s] (4 - s) 


Moreover, it can be checked that these solutions are attractors of neighbouring cosmological 
flows. As a result, the brane evolution is always given by the above equations after a sufficiently 
long time, independently of the initial conditions. A characteristic property of the solutions is 
that C 2 < 0 and, therefore, x < 0. On the other hand, Ci + C 2 > 0 and the Hubble parameter 
is always real. 

The solutions (21) describe an expanding Universe with an acceleration parameter 


Qo 



(25) 


This expression is independent of the equation of state of the brane matter (as long as tc 7 ^ 1/3). 
The Hubble parameter, however, depends on w, as can be seen from eqs. (16), (24). For q < 3 
the parameter Qo is positive and the expansion is accelerating. 

For 0 < g < 4 the scale factor increases as a power of t for long times 


do 


+ |(C^l + ^^2) 


1/2 


{t - ti) 


3/q 


(26) 


For q = 0 the expansion is exponential 


do 


- U/02C- , 


(27) 


independently of the equation of state of the brane matter. This case corresponds to the hxed- 
point solutions studied in a more general context in ref. [ 8 ]. It is remarkable that one can obtain 
exponential expansion even for non-relativistic brane matter. However, the value q = 0 implies 
the presence of vacuum energy density in the bulk, which is assumed to be transferred to the 
brane particles. 

The case of relativistic brane matter must be considered separately, as eqs. (16), (21), (24) 
give Ho 0 for w 1/3. The solution (21)-(24) is approached only after a very long time in 
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this limit and the initial conditions determine most of the evolution. For w = l/2> the solution 
is 


1 

6 MI 


p{ao) 


x{ao) 


1 Pjati , T 1 2-q 

6M|,; 6-q 

Xiati _ T 1 2-q 

at ’ 


(28) 

(29) 


where aoi = ao{ti) etc. The sum p/{6Mpi) + x that determines the Hubble parameter behaves 
as an effective radiation term ~ a“^. This is achieved through the cancellation of the dominant 
second terms in the r.h.s. of eqs. (28) and (29). As a result, the Universe is always decelerating 
for w = 1/3. 


Discussion: Our main result is summarized by eq. (25); The influx of energy onto the 3- 
brane from the bulk can lead to accelerated expansion, depending on the equation of state of 
the bulk matter. For positive acceleration one needs q < 3 for the variable q that parametrizes 
the dependence of the energy flow on the scale factor (see eq. (20)). If the bulk is populated 
by a gas of non-relativistic particles that drift slowly towards the brane, we expect q = 3. If 
the drift velocity increases with time, q < 3. A value of q slightly below 3 leads to a positive 
acceleration parameter smaller than 1, the case favoured by the astronomical data. 

The construction of a complete model that realizes this scenario requires the technically 
difficult solution of the Einstein equations in the bulk. However, certain properties of such a 
solution can be inferred in general terms. The energy flow along the fifth dimension towards 
the brane implies the presence of a region of high energy density where the flow originates. For 
example, one could consider a second brane that emits energy into the bulk. Its evolution is 
governed by eqs. (15)”(17) with a positive A solution analogous to that of eqs. (21)” 
(24) exists, with H < 0. It predicts a Big Crunch after a finite time. However, neighbouring 
cosmological flows diverge from this solution. There are also various solutions with H > 0. The 
generic cosmological evolution leads to the depletion of the energy density on the second brane. 
The flow of energy along the fifth dimension is expected to stop after a certain period, whose 
length depends on the initial energy density of the brane. 

These examples indicate that, in a complete model, the solutions (21)-(24) are expected to 
be valid only for a finite time interval. This makes the derivation of exact analytical expressions 
more difficult. However, we believe that the physical arguments for the form of the energy- 
momentum tensor that we gave in the previous section (see eq. (20)) give a valid approximation, 
especially in the case of slow accretion of energy by the brane. 

The RS model is employed here only in order to provide the context in which our arguments 
can be implemented. The discussion of cosmological solutions must take into account the grav¬ 
itational effects of the brane tension (the vacuum energy associated with the brane) and the 
matter localized on it. The RS model is the only example in which this has been realized in such 
a way that conventional cosmology is reproduced in the low density limit. However, we view 
this model only as a toy one. Many of its ingredients, such as the negative energy density in 
the bulk and the structure of the resulting AdS space, do not seem crucial for the emergence of 
accelerating solutions. This view is supported by the discussion of eq. (1) in the introduction. 
In order to test this conclusion, it would be interesting to implement this mechanism in a model 
with a bulk geometry that is flat, at least asymptotically. 
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An important property of the realization of the mechanism in the context of the RS model 
is the absence of acceleration during radiation domination, even for energy influx. For brane 
matter with re = 1/3 the solution (28), (29) describes a radiation dominated brane, whose 
expansion is decelerating. This means that the cosmic acceleration is a phenomenon only of the 
era of matter domination. 

It must be pointed out that the solution (28), (29) depends crucially on the presence of the 
component y of eq. (12). In the case of the pure AdS bulk space of the RS model this component 
describes, according to the AdS/CFT correspondence, the boundary conformal field theory [11]. 
It has been characterized as “mirage” or “Weyl radiation”. In our scenario we have assumed a 
more general matter content for the bulk theory, for which this interpretation is not established. 
However, the term ~ ^Hx in the evolution equation of y is characteristic of radiation. It is 
not clear if a similar “radiation” component will be present in a different realization, with an 
asymptotically flat bulk space for example. Therefore, the absence of acceleration during the 
radiation dominated era on the brane, even for energy influx, may be a property of RS-type 
models only. 

The energy influx may not be present at all times for other reasons. For example, it may be 
significant only for very low energy densities on the brane, and be replaced by energy outflow for 
densities exceeding a critical value. This again would lead to the conclusion that the accelerated 
expansion is only a recent phenomenon in cosmic terms. The dynamical localization mechanism, 
that determines the details of energy accretion, is the main direction of further research in order 
to understand this point and construct phenomenologically viable models. 


References 

[1] J.L. Sievers et at, arXiv:astro-ph/0205387; 

A.G. Riess et at, Astron. J. 116 (1998) 1009 [arXiv;astro-ph/9805201]; 

S. Perlmutter et at, Astrophys. J. 517 (1999) 565 [arXiv:astro-ph/9812133]; 

W.J. Percival et at, Mon. Not. Roy. Astron. Soc. 327 (2001) 1297 [arXiv:astro- 
ph/0105252]. 

[2] For a review, see B. Ratra and P.J. Peebles, arXiv:astro-ph/0207347. 

[3] F. Hoyle, Mon. Not. Roy. Astron. Soc. 108 (1948) 372; 109 (1949) 365; 

F. Hoyle and J.V. Narlikar, Proc. Royal Soc. A 270 (1962) 334; 277 (1964) 1; Ann. Phys. 
54 (1969) 207. 

[4] L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 3370 [arXiv:hep-th/9905221]; 
Phys. Rev. Lett. 83 (1999) 4690 [arXiv;hep-th/9906064]. 

[5] C. Csaki, M. Graesser, G.F. Kolda and J. Terning, Phys. Lett. B 462 (1999) 34 [arXiv:hep- 
ph/9906513]; 

J.M. Gline, G. Grojean and G. Servant, Phys. Rev. Lett. 83 (1999) 4245 [arXiv:hep- 
ph/9906523]. 

[6] P. Binetruy, G. Deffayet and D. Langlois, Nucl. Phys. B 565 (2000) 269 [arXiv:hep- 
th/9905012]; 

P. Binetruy, C. Deffayet, U. Ellwanger and D. Langlois, Phys. Lett. B 477 (2000) 285 
[arXiv;hep-th/9910219]. 

[7] C. van de Bruck, M. Dorca, C.J. Martins and M. Parry, Phys. Lett. B 495 (2000) 183 
[arXiv;hep-th /0009056]. 


7 



[8] E. Kiritsis, G. Kofinas, N. Tetradis, T. Tomaras and V. Zarikas, arXiv;hep-th/0207060. 

[9] A. Hebecker and J. March-Russel, Nucl. Phys. B 608 (2001) 375 [arXiv;hep-ph/0103214]; 
D. Langlois, L. Sorbo and M. Rodriguez-Martinez Phys. Rev. Lett. 89 (2002) 171301 
[arXiv;hep-th /0206146]. 

[10] R. Jackiw and C. Rebbi, Phys. Rev. D 13, 3398 (1976); E. Weinberg, Phys. Rev. D 24, 
2669 (1981); V.A. Rubakov and M.E. Shaposhnikov, Phys. Lett. B 125, 136 (1983); A. 
Barnaveli and O. Kancheli, Sov. J. Nucl. Phys.51, 573 (1990); ibid. 52, 576 (1990); G. 
Dvali and M. Shifman, Phys. Lett. B 396, 64 (1997). 

[11] S.S. Gubser, Phys. Rev. D 63 (2001) 084017 [arXiv;hep-th/9912001]. 



